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The purpose of this paper is to study the dynamics of a square bilHard with a non-standard reflection law such 
that the angle of reflection of the particle is a linear contraction of the angle of incidence. We present numerical 
and analytical arguments that the nonwandering set of this billiard decomposes into three invariant sets, a 
parabolic attractor, a chaotic attractor and a set consisting of several horseshoes. This scenario implies the 
positivity of the topological entropy of the billiard, a property that is in sharp contrast with the integrability 
of the square billiard with the standard reflection law. 
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A billiard is a mechanical system consisting of a 
point-particle moving freely inside a planar region 
and being reflected off the perimeter of the region 
according to some reflection law. The specular 
reflection law is the familiar rule that prescribes 
the equality of the angles of incidence and reflec- 
tion. Billiards with this reflection law are conser- 
vative systems, and as such are models for phys- 
ical systems with elastic collisions. For this rea- 
son and their intrinsic mathematical interest, con- 
servative billiards have been extensively studied. 
Much less studied are dissipative billiards, which 
originate from reflection laws requiring that the 
angle of reflection is a contraction of the angle of 
incidence. These billiards do not preserve the Li- 
ouville measure, and therefore can model physical 
systems with non-elastic collisions. In this pa- 
per, we investigate numerically and analytically 
a dissipative billiard in a square. We find that 
its dynamics differs strikingly from the one of its 
conservative counterpart, which is well known to 
be integrable. Indeed, our results show that a dis- 
sipative billiard in a square has a rich dynamics 
with horseshoes and attractors of parabolic and 
hyperbolic type coexisting simultaneously. 
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I. INTRODUCTION 

Billiards are among the most studied dynamical sys- 
tems for two main reasons. Firstly, they serve as models 
for important physical systems (see e.g. the book^ and 
references therein), and secondly, despite their simplic- 
ity, they can display a rich variety of dynamics ranging 
from integrability to complete chaoticity. Most of the ex- 
isting literature on billiards is devoted to billiards with 
the standard reflection law: the angle of reflection of the 
particle equals the angle of incidence (cf.''^-^*'). These bil- 
liards are conservative systems, i.e. they admit an invari- 
ant measure that is absolutely continuous with respect to 
the phase space volume. 

In this paper, we are concerned with billiards with a 
non-standard reflection law according to which the angle 
of reflection equals the angle of incidence times a constant 
factor < A < 1. Since we have observed numerically 
that such a law has the effect of contracting the phase 
space volume, billiards with this law will be called "dis- 
sipative" in this paper. 

Recently, Markarian, Pujals and Sambarino^ proved 
that dissipative planar billiards (called "pinball billiards" 
in their paper) have two invariant directions such that the 
growth rate along one direction dominates uniformly the 
growth rate along the other direction. This property is 
called dominated splitting, and is weaker than hyperbol- 
icity, which requires one growth rate to be greater than 
one, and the other one to be smaller than one. The result 
of Markarian, Pujals and Sambarino applies to billiards 
in regions of different shapes. In particular, it applies to 
billiards in polygons. This is an interesting fact because 
the dominated splitting property enjoyed by the dissi- 
pative polygonal billiards contrasts with the parabolic 
dynamics observed in the conservative case^iS. 

Here, we take a further step towards the study of 
dissipative polygonal billiards analyzing the dissipative 
square billiard. Taking into account the symmetries of 
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the square, we perform our analysis on a reduced phase 
space. We provide theoretical arguments and numerical 
evidence that the nonwandering set of our system de- 
composes into three (possibly empty) invariant sets: a 
parabolic attractor, a hyperbolic attractor and a horse- 
shoe. This dynamics is clearly richer than the one of the 
conservative square billiard, which is a fully integrable 
system. In this paper, we also conduct a rather detailed 
numerical study of the changes in the properties of the 
nonwandering set as the parameter A varies. 

We should mention that results somewhat similar to 
ours were obtained for non-polygonal billiardsi^ and the 
dissipative equilateral triangle billiard^. 

The paper is organized as follows. In Section |lll we 
give a detailed description of the map for the dissipative 
square billiard. Some results concerning the invariant 
sets of this map are presented in Section Hill To study our 
map, it is convenient to quotient it by the symmetries of 
the square. This procedure is described in Section IIVI 
and produces the so-called reduced billiard map. Sec- 
tion |V] is devoted to the study of two families of periodic 
points of the reduced billiard map. In particular, we show 
the stable and unstable manifolds of a fixed point of the 
reduced billiard map (corresponding to a special periodic 
orbit of the billiard map) have transversal homoclinic in- 
tersections, and use this fact to conclude that the dis- 
sipative square billiard has positive topological entropy. 
Finally, Section IVII contains the bifurcation analysis of 
the nonwandering set of the reduced billiard map. 



II. THE SQUARE BILLIARD 

Consider the square D — [0, 1] x [0, 1] C M^. For our 
purposes, D is called the square billiard table. To study 
the dynamics of the billiard inside this table, it is suffi- 
cient to know the angle of incidence at the impact points 
and the reflection law. For the usual reflection law (the 
angle of reflection is equal to the angle of incidence) we 
find the next impact point s' and angle of reflection d' 
by the billiard map (s', B') — B{s, 9) acting on the previ- 
ous impact {s,0). This map admits an explicit analytic 
description. Its domain coincides with the rectangle 



M ^ [0,4] X 
from which the set 
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S+ = { {s,d) e X: {s} = or {s} + tan6l e {0,1}} 

is removed. The symbols [s] and {s} — s — [s] stand 
for the integer part and the fractional part of s, respec- 
tively. An element of corresponds to an orbit leaving 
or reaching a corner of D (see Fig. [1]) . By reversing the 
role of time in this description of S'^ , one obtains the set 

S- = {{s,9) e M: {s} = or (s, -X-'^O) £ S+ } . 




FIG. 1: Phase space M\S^ 



Both sets and S consist of finitely many analytic 
curves. Next, let 



Ml - {(5,6*) eM: 
M2 = {(3,9) eM: 



I: {s} > and {s} -I- tan 6* > 1 } , 
! : {s} > and < {s} + tan6l < 1 } , 
! : {s} > and {s} + tan6l < } . 

The bihiard map B: M\S+ ^ M\S- is defined by 

1-H 



B(s, 



's] + 1 



9 on Ml, 



(W - 1 



tan^ ("^°^^)'2 
{s} - tan 6* (mod 4), -6*) on A^2, 

on M3. 



(mod4),-^- 
tan^ ^ ' 2 



This map is clearly an analytic diffeomorphisni in its do- 
main. The inverse of B is easily obtained by noticing 
that the billiard map is time-reversible. That is, given 
the map T(s, 9) — (s, —9), we have 

B-^ =ToBoT-\ 

To modify the reflection law, we compose B with an- 
other map TZ: M ^ Ai. The resulting map ^ = TZ o B 
is called a billiard map with a modified reflection law. 

Several reflections laws have been consideredi*2i. In 
this paper, we consider the following "dissipative" law. 
Given < A < 1, we set 

nx{s,9) = is,X9). 

According to this law, the direction of motion of the par- 
ticle after a reflection gets closer to the normal of the 
perimeter of the square (see Fig. [2]). To emphasize the 
dependence of the billiard map on the parameter A, we 
write 

$A = 7^A o B. 

As a side remark, one can also define the map $a for 
A > 1. In this case, the map TZ\ expands uniformly the 
angle 9, and $a becomes a map with holes in the phase 
space. It is interesting to observe that the maps 
and are conjugated for < A < 1. Indeed, it is not 
difficult to check that 

$A-i = (7^A o T)-' o $-1 o (7^A o r), 
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FIG. 2: Dissipative reflection law. 



by using the fact that T and IZx commute and that 
TL'^ = Tix-i- Therefore, all the results presented in 
this paper hold for A > 1 as well, provided that we re- 
place the word "attractor" with the word "repeller" , and 
switch the words "stable" and "unstable" . 



III. HYPERBOLICITY 

Let {so,0o) be an element of M \ S^. Set (si,0i) ~ 
^x{so, do), and denote by t{so, Oq) the length of the seg- 
ment connecting sq and Si. Using elementary trigonom- 
etry, one can show in a straightforward manner that the 
derivative of $a takes the following form: 



cos 6>o t{so,0o) 
cosA~^6'i cosX^^9 







A 



In fact, the previous formula holds for every polygon, and 
not just for the square (see Formula 2.26 in ref.^). 

Now, suppose that {(si,0i)}"^Q are n + 1 consecutive 
iterates of Then, we see that 



i?$^(so,0o) = (-l) 
where 

cos 9o 



Pn{so,Oo) 



and 



n 



cos 9i 
cos X^^9i 



Mso,Oo) = X" 



Cn(so, 



1 



n-1 



cos A ^9. 



k=i+l 



COS 9k 
cos X^^9k 



We now prove a simple lemma concerning the stability 
of the periodic points of $a- It is not difficult to see 
that this result remains valid for every polygon and for 
other reflection laws (e.g. TZ\{s,d) = {s,9 — csm29) with 
< c < 1/2 as in^). 

Lemma III.l. For every X e (0, 1), the periodic points 
of^x of period 2 and period greater than 2 are parabolic 
and hyperbolic, respectively. 



Proof. Suppose that (sq, ^o) is a periodic point of <I>a with 
period n. Since (s„,0„) = {so,9q), it turns out that 



an{so,9o) = 



i=0 



cos 9i 
cos A^^ 



Now, note that each term cos 9i/ cos X~^9i in the ex- 
pression of Q!„(soj^o) is equal or greater than 1 with 
equality if and only if 9i — 0. Also, note that D^^{so,9o) 
is a triangular matrix, and so the moduli of its eigenval- 
ues are a„(so,6'o) and A" < 1. Therefore to determine 
the stability of (so,0o) is enough to find out whether or 
not an{so,9o) is greater than 1. 

If n = 2, it is easy to see that the trajectory of 
{sq,9o) must always hit the boundary of D perpendic- 
ularly. In other words, we have 9q = 9i = 92 = 0, and so 
^2 (30,6*0) = 1- Periodic points of period 2 are therefore 
parabolic. Clearly, a necessary condition for a polygon 
to admit periodic points of period 2 is that the polygon 
must have at least 2 parallel sides (not a sufficient con- 
dition though). 

Now, suppose that n > 2. In this case, we claim 
that (soj^o) is hyperbolic. Indeed, when n > 2, the bil- 
liard trajectory of (so,0o) must have at least two non- 
perpendicular collisions with the boundary of 13, and 
since cos 9i/ cos X^^9i > 1 for such collisions, we can im- 
mediately conclude that q;„(so,^o) > 1- CH 

A more elaborated analysis along the lines of the proof 
of Lemma IIII.ll yields some general conclusions on the 
chaotic behavior of general dissipative polygonal bil- 
liards. Two of such conclusions are stated in Proposi- 
tions nil. 21 and IIII.3I below. To state these proposition, 
we need first to introduce the notion of uniformly and 
non-uniformly hyperbolic systems. Unless specified oth- 
erwise, $A denotes the map of a dissipative billiard in a 
general polygon D throughout the rest of this section. 

A set S C is called invariant if <I>^^(S) = S. An 
invariant set E is called hyperbolic if there exist a norm 
II ■ II on A^, a non-trivial invariant measurable splitting 
T^M — E^®E^ and two measurable functions < /i < 1 
and K > Q ox\Yj such that for every (s, ^) G E and every 
n>l, we have 

\mi\E^f^s,o)\\<K{s,9)ii{s,er, 

If the functions /i and K can be replaced by constants, 
then E is called uniformly hyperbolic, otherwise it is 
called non-uniformly hyperbolic. 

We can now state our propositions. The first one con- 
cerns billiards in polygons without parallel sides. For 
such polygons, the map $a does not have periodic points 
of period 2. 

Proposition III. 2. Let D be a polygon without parallel 
sides, and suppose that E is an invariant set of^x- Then 
E is uniformly hyperbolic for every X G (0, 1). 
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The second proposition concerns billiards in rectangles. 
In this case, has periodic points of period 2. Denote 
by V the set of all these points. It is not difficult to check 
that 7-" is a parabolic attractor for every A G (0, 1). 

Proposition III. 3. Let D he a rectangle, and suppose 
that S is an invariant set of^x not intersecting the basin 
of attraction ofV. Then there exists A, £ (0, 1) such that 
E is hyperbolic for every A G (0,A,), and is uniformly 
hyperbolic for every A G (A*, 1). 

For reasons of space, the proofs of these propositions 
are omitted and will appear elsewhere^. 



IV. THE REDUCED BILLIARD MAP 




FIG. 3: Periodic orbit. 



The analysis of the billiard dynamics can be simplified 
if we reduce the phase space. First, we identify all sides 
of the square by taking the quotient with the translations 
by integers of the s-component. Then, due to the symme- 
try along the vertical axis at the midpoint of the square, 
we can also identify each point (s, 6) with (1 — s, —9). To 
formulate the reducing procedure more precisely, we de- 
fine an equivalence relation ~ on by (si, ^i) (s2, ^2) 
if and only if 7r(si,6'i) = 7r(s2,^2): where tt: — > is 
the function defined by 



Ti(s,e) 



if e G [0, f ) , 
if^G (-f,o) 



Let M denote the image of tt. Clearly, we have 



M = (0, 1) X [0, 1) 



Note that it is possible to identify the set M with the 
quotient space A4/ ~. We call M the reduced phase 
space. The induced billiard map on M is the reduced 
map, which we will denote by (px. 

It is clear from the definition of tt that tt^^{s,9) con- 
sists of 8 elements for every {s,9) G Af, and so (A^,7r) 
is an 8-fold covering of M . It is then easy to see that 
the reduced billiard map (px is a factor of the original 
billiard map $a by noting that the quotient map tt is 
indeed a semiconjugacy between (f>x and $a, i-e. we have 

that TT O $_)^ = O TT. 

In what concerns the relation between the dynami- 
cal systems defined by $a and there are several key 
points that are worth remarking. First, we note that pe- 
riodic points of 4>x lift to periodic points of $a- To be 
more precise, an orbit of period n under (px is lifted to 
either eight orbits of period n, or four orbits of period 2n, 
or two orbits of period 4n, or one orbit of period 8ri for 
$A- Analogous statements can be produced for the lifts 
of transitive sets and the existence of invariant measures. 
Namely, transitive sets for cjix are lifted to a finite num- 
ber of transitive sets for and any invariant measure 
under the dynamics of 4>x corresponds to a finite number 
of invariant measures under . Finally, we remark that 



the reduced map <j)x has positive topological entropy if 
and only if this is the case for the billiard map $a- 

By studying the trajectories of the billiard map we 
have basically two cases: either the billiard orbit hits 
a neighboring side of the square or the opposite side. 
Separating these cases there is a corner which is reachable 
only if the initial position (s, 9) G M is in the singular 
curve 

S+ = {{s,9) G M: s + tan^^ 1}. 

This curve separates the reduced phase space in two con- 
nected open sets: Mi below 5*+ and M2 above 5*+. 

Let /i : Ml M and /2 : M2 -> M be the transfor- 
mations defined by 

/i(s, 9) = {s + tan0, \9) for (s, 9) G Mi, 

f2{s,9) = ((l-s)cot0,A(|-e)) for {s,9) ^ A'h. 

The reduced billiard map is then given by 

on Ml, 
^ [/2 onM2. 

Its domain and range are M\S^ and M\S' 
where 



respectively. 



{{s,9) ^ M: s- tan(A" 



0} 



Like the billiard map $a, the reduced billiard map 0a is 
an analytic diffeomorphism. Notice that (px maps hor- 
izontal lines into horizontal lines, a consequence of the 
fact that its second component is independent of s. 
Finally, we observe that the subsets of M where the 



maps 
spectively. 



and 



are defined for every n > are, re- 



M+ = M \ y (j}l'^{S+) and 



M- 



n>0 
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V. ATTRACTORS AND HORSESHOES 

We start this section by formulating several definitions. 
The stable set of an element g G A/ is defined by 



Wiq) = \ xeM 



liin ||01H^)-0A 







where || • || is the Euclidean norm on M. In the case of 
an invariant set A = (j)\{A), we define its stable set to be 

The unstable sets W'^{q) and VF"(A) are defined anal- 
ogously by replacing 0a with and M+ with M^. 
When W'^^'^^A) turns out to be a manifold, we will call 
it an unstable{stable) manifold. 

Suppose that A is an invariant subset of M . Then A 
is called an attractor if A = W'^{A) and WiA) is open 
in M+, and is called a horseshoe if neither W'(A) is an 
open set in M+ nor is W"{A) an open set in M~. Note 
that a saddle periodic orbit is a horseshoe according to 
this definition. A finite union of hyperbolic invariant sets 
Ai, . . . , Am is called a hyperbolic chain if 

W{A,)r]W''{A,+i) / for i = l,...,m- 1. 

A point X £ M-|_ is said to be nonwandering if for ev- 
ery open neighborhood U containing there exists n>l 
such that C/n0^(C/) 7^ 0. We denote by C Af + flAf" 
the set of all nonwandering points of (j)x ■ We say that two 
hyperbolic periodic points x, y G ^x are homoclinically 
related if W^{x) and W'^{y) intersect transversally, and 
W^{y) and W^{x) intersect transversally. The closure in 
Af + n of the set of periodic points homoclinically re- 
lated to a hyperbolic periodic point x G J^a is called the 
homoclinic class of x. Every homoclinic class is a transi- 
tive invariant subset oiftx (see Refi^ (Ch. IX, Prop. 5.2)). 



A. Parabolic attractor 

Let us define P = t^{V). It is easy to see that 

P = {{s,e) G M: e' = o}, 

and each point of P is a parabolic fixed point coming 
from period 2 orbits of the original billiard (orbits that 
bounce between parallel sides of the square). It is an 
attractor and W^*(P) includes the set of points B that 
are below the forward invariant curve 

Soo = |(s,6') e M: s-^-^tan(A*6l) = l| . 

The sequence (t)1{Soo) converges to the point (1,0). The 
pre-image of B is at the top of phase space. Moreover, 
its basin of attraction is 



TJ/'S 


















B ^ 







FIG. 4: Invariant manifolds of and singular curves 
for the reduced billiard map (A = 0.6218) 



By Proposition IIII.3| the set ^lx\ P 'is hyperbohc. 
Hence, every periodic point in JIa \ P has stable and 
unstable manifolds. Because of the cutting effect of the 
singular sets S~ and S'^ , these manifolds are countable 
unions of smooth curves. 



B. Fixed point and its invariant manifolds 

The map $a has many periodic orbits. Two special 
periodic orbits of period 4 can be found by using the 
following simple argument. A simple computation shows 
that if an orbit hits two adjacent sides of the square with 
the same reflection angle 6'a, then 

ttA 



2(1 + A) 



If we further impose the condition that the orbit hits the 
two sides at si and S2 in such a way that {si} = {52} — 
sx, then we obtain 



1 



Sx 



1 + tan Ox 



n>0 



By symmetry, we conclude that ^\{sx,Ox) ~ {sx,0x)- 
Using once again the symmetry of the square, we also 
have $^(1 — sa, — ^a) = (1 — sa,— 6'a)- One of these orbits 
is depicted in Fig. |3l 

Due to the phase space reduction, the periodic orbits 
just described correspond to the fixed point 

Pa = isx,Ox) 

of (f>x ■ This is actually the only fixed point of <j)x in M2 
outside of P. By Lemma rill.li px is hyperbolic and thus 
it has local stable and unstable manifolds W^^^{px) for 
every A G (0,1). Since 4>x maps horizontal lines into 
horizontal lines, and the set S~ does not intersect the 
horizontal line through px , we see that the local unstable 
manifold of px is given by 

W,lM = {{s,e)eM:9^9x}. 
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In fact, the global unstable manifold consists of a collec- 
tion of horizontal lines cut by the images of S*^. 

The geometry of the stable manifold is more compli- 
cated. By definition points on the stable manifold con- 
verge to the fixed point. Moreover, Wf^^{p\) cannot cross 
S^. Thus, W{^^{p\) is contained in M2. The graph trans- 
form associated with the corresponding branch of (j)x is 
the transformation 

r{h){e) = 1 - h{gx{e)) tan^, 

where g\: [0,7r/2) — > [0,7r/2) denotes the affine contrac- 
tion 

Iterating k times the zero function by T we obtain 

k-l ji-l 

Hence, the local stable manifold of px is the curve 
WiUpx) = { {hx{0), 9):0<e<^and0<hx{9)<l}, 
where 

00 n — 1 

hxi9)^Y.(-^)"]l'M9m)- (1) 

n=0 i=0 

This series converges uniformly and absolutely since 
tan((;^'(6')) converges to tanOx as n — > 00, and < 
tan6'A < 1. The same statement holds for the series of 
the derivatives of hx- Thus, hx is smooth. 

The invariant manifolds of px, the singular curves of 
the reduced billiard map and the upper boundary of 
B are depicted in Fig. 2] 

Let A2 be the unique solution of 

/ia(A6Ia) ==tan(6lA) for A G (0, 1) . 

Geometrically, A2 is the value of A such that the singu- 
lar set S~, the local manifold Wf^^{px) and the closure 
of the first iterate of the piece of the unstable manifold 
of Pa contained in Mi have non-empty intersection (see 
Fig. [3b)). A numerical computation shows that 

A2 = 0.8736... 

Let A be the closed set bounded by Wf^^{px) and 
W^IM as in Fig. [Sib). 

Proposition V.l. For A > A2, there is a compact (fix- 
invariant set Aq C int(A) such that flxHA C {px}U Aq. 

Proof. For A > A2, we have /i(AnMi) C int(A), and wc 
can a find a compact forward- invariant set Aq C int(A) 
under /2 such that /i(A n Mi) C Aq. Hence Aq is also 
(/)A -invariant. Since px is the only nonwandering point 
in A whose orbit does enter Mi, we have J7a H A C 
{pa}UAo. □ 



Proposition V.2. The invariant manifolds of px have 
transverse homoclinic points if and only i/ < A < A2 . 

Proof. To prove the existence of homoclinic points we it- 
erate a piece of the local unstable manifold in Mi and 
show that it intersects transversely the local stable man- 
ifold in M2. Taking into account that 0a maps horizontal 
lines into horizontal lines, 4>x\mi ~ /i and 

/i(iyi^, n Ml) = (tan^A, 1) X {XOx} 

the problem of finding homoclinic intersections reduces 
to proving the following chain of inequalities: 

tan6iA < hx{\Bx) < 1. 

We will see that these inequalities hold if and only if 
< A < A2. 

Lemma V.3. The inequality tan9x < hx{X9x) holds if 
and only i/ < A < A2 . 

Proof. For A > sufficiently small we have tan 0a + 
tan(A6'A) < 1 since 6'a as A — ?> 0. On the other hand, 
we know by definition of hx that hx{X9x) + tan(A6'A) > 1 
for every A S (0, 1). Putting these two inequalities to- 
gether we conclude that 

tan 61a < 1 - tan(A6lA) < hxiX9x) 

for every A > sufficiently small. Since hxiX9x) — tan^A 
is strictly decreasing for A € (0, 1) and A2 is the unique 
solution of hx{X9x) = tan6'A we obtain the desired result. 

□ 

Lemma V.4. The inequality hx{X9x) < 1 holds for ev- 
ery X e (0, 1). 

Proof. At the fixed point we compute 

^''^^''^ ^ "(l-AtanTA)a + tan0A) ^ ' 
Define now 

m{9) = hx{9x) + h'x{9x){9-9x) 

_ 1 - A tan 61a - sec^ 6Ia [9 - 9x) 
~ (1 - A tan6'A) {l + tim9x) 

to be the function whose graph is the line tangent to the 
graph of hx at the fixed point px- Since hx is concave, 
hx{9) < m{9) for every 9 G (0,7r/2). Thus, it is enough 
to check that m{X9x) < 1 for A < 1, by using elementary 
estimates. □ 

Thus, fiiW^lM n Ml) n Wf^M ^ if and only if 
< A < A2. To conclude the proof of the proposition, 
note that for A > A2 , by Proposition IV. II the region A 
is a trapping set. Since W^{px) C int(A) U W^^{px), px 
has no homoclinic intersections for A > A2. □ 

The following corollary is a direct consequence of 
Proposition IVj and Ref.-^ (Ch. 7,Th. 4.5). 

Corollary V.5. The map (px has positive topological en- 
tropy for every < A < A2 . 
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(a) A = Ai . Since the image of the map (^-\ is always below the 
line 6 = A7r/2, the basin of attraction of P is only B U cj)~^^{B). 
Moreover, the region in light gray is forward-invariant. 



(b) A = A2. The shaded region A between the stable and unstable 
local manifolds of p\ is forward-invariant. 



FIG. 5: Trapping regions for 



C. Two Families of Periodic Orbits 

Given (n,m) G N^, a straightforward computation 
shows that 

° /r(s, 0) = ((-i)"-iT„,„(0, s), gi{y^e)) , 

where T„^m is given by 

and the sequences of functions hn , 7n and 5„ are defined 

by 

n i— 1 

in{e) = n ^oi{g\{e)) 

i=0 



5„(e') =^tan(A* 



1=0 

Recall that g\ is the affine contraction g\{9) — \{tt/2~9). 
For each n>\, define g„ and Pn as the unique solutions, 
when they exist, of 

fl ° fi{(ln) = qn and /f'-^ o /i(p„) = p„ . 

In agreement with this definition, we set go = Pa- 

Proposition V.6. There exists a unique decreasing se- 
quence Cn G (0,A2) such that qn is a periodic point of 
period n + 2 for ipx if and only if X G (0, c„). 



Proof. Let qn = (s„,0„). A simple computation shows 
that 



s„ - (1 - SM) 

n ^A(l-A) 



72(A"g„) 
72(A"0. 



1' 



(3) 



" 2 1 - A2+" ■ 

Let A„ be the set of points (s. G Afi such that 

l-^„(0) < S < 1 ~ Sn-i{e) . 

Since /"(A„) C A/2, to show that exists it is enough 
to check that {sn,dn) belongs to A„. 

In the following estimates we will frequently use the 
fact that the tangent is a convex function, i.e., for every 
x,y & (0, ^) and < A < 1 we have 

tan(Aa; + (1 — X)y) < Atan(.T) + (1 — A)tan(j/) , 

and moreover that for every < x < tan(a;) > x. 

Let us start by proving that s„ > 1 — Sn{0n)- By 
definition 

72(A"0„) = cot(A"0„) cot (a (I - A"0„)) 
cot(A"6i„) 



> 



Atan(f - A'' 



Thus 72 (A"^^) > 1 which, taking into account the def- 
inition of s„, implies that s„ > 1 — S'„(0„). To prove the 
other inequality, we start by noting that 



l-5„-i( 



l-^n-l(^n)-COt(gA(A"g„)) 
72(A"0„)-1 



(a) qi 



(b) q-2 




(d) pi (c) P2 (f) P3 

FIG. 6: Billiard trajectories corresponding to the periodic points qn and p„ for A = 0.6. 



Since 72 (A" (^„) > 1, we only need to prove that 

5„_i(0„)+cot(5A(A"e„)) > 1. (4) 
Using the definition of 5'„ we get 



"~ . 1 - A 

Sn—l 

(^n) = ^tan(y^„)>^„- 

i=0 



A 



On the other hand 



cot(5A(A"^„)) = tan(A-ie„) > -j- . 

Putting these estimates together we obtain 

1 - A"+i 



^n-l(^n)+COt(^?A(A"^^„)) 



A(l - A) 



TT 1 - A"+^ TT n + 1 

2 1 - A"+2 ^ 2 n + 2 



> 1 



for every A e (0, 1). 

It remains to prove that s„ > 0. It is clear 
that dSn{9n{\)) / d\ > for every A e (0,1). Since 
Sn{Onm = and 5„(0„(1)) = (n + 1) tan(7r(n + 
2)^72) > 1 for n > 1, wo conclude that S'„(6'„(A)) < 1 
if and only if A G (0,c„). Here c„ G (0, 1) is the unique 
solution of 



^tan(y^„(A)) =1. 



i=0 



Thus s„ > 1 - S„(6»„(A)) > if and only if A G 

(0,c„). Now wc prove that {c„} is decreasing. Since 
dSn{dnW) / dX > it is sufficient to prove that 

) > Sn{0n) . 



By definition of 5„ we have that 

n 

= Sn{en) + J2 [tan(V^„+i) - tan(A'e„)] 

+ tan(A"+i(?„+i). 

Let AS'„ = 5„+i(6'„+i) - 5„(6i„). Note that 6i„ > 6i„+i. 
Since tan(y — a;) > (tan(2/) — tan(a;))(l — tan(a;) tan(y)) 
for every < a; < y < 7r/4 we get 

A<? ^tflna^+ifl ^ ^ tan(A'(6>„ -6i„+i)) 
A5„>tan(A " X. 1 _ tan(A^^„+:) tan(A^^„) 



i=0 



-><-PT,a"+ifl ^ >^ tan(A'(6'„ - On+i)) 
>tan(A ^n+i)-l^ i_tan^(^„) " 

In the derivation of the previous inequality we have used 
the upperbound: tan(A*0„4-i) tan(A*0„) < tan^(0„) for 
every n > 1. Let 



Pn 



A"+2 



1 + A + • • ■ + A"+i ■ 
Clearly 6n = (1 + Since p„/A"+^ < 1 we obtain 



A5„>tan(A"+i0„+i)-^^ 

i=0 



tan( p„A'6i„+i) 
•tan2(6i„) 



>tan(A ^n+i) ^l-An+i(l-tan^(M(l-A) 



Using the expression for p„ we get 

ASr, 



^ tan(A"+i0„+i)(l - A) / _ ^Atan2(^\ 
'"^ (l-tan2(0„))(l-A"+2) 20n J 
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Since < 9n < Oi < 7r/4 we have that 
tan2(0„) < (^^^ tan2((?i) 

< ^(1 + A + A^) tan^ ( ^^^^ \ A 
ttA ^ ^ V2A(1 + A + A2)y 

Hence A^n > as we wanted to show. Finally, {c„} is 
bounded from above by A2 since a numerical computation 
reveals that 

ci = 0.7964... 

which, taking into account the numerical value of A2, 
implies that c„ < ci < A2. □ 

The proof of the next result is omitted because it is 
similar to the previous one. 

Proposition V.7. // A £ (0, A2) then pn is a periodic 
point of period 2n for 4>x . 

By Lemma IlII. 11 these periodic points are hyperbolic. 
As we shall see in the next section, these orbits seem 
to play an important role in the dynamics of 4>\ for dif- 
ferent values of A. The corresponding billiard orbits in 
configuration space are depicted in Fig. |6l 

Since the sequence c„ is decreasing, we can define Ai = 
lim„^ooC„. The number Ai is also the unique solution 
of the equation 

00 

^ tan (|a'+1(1 - A)) = 1 for A e (0, 1). 

In geometrical terms, when A = Ai, the intersection point 
of the curve S^o with the line 9 = 7rA(l — A)/2 lies ex- 
actly on the line s = (see Fig. [3a))- This intersection 
point is also the limit of the sequence g^. A numerical 
computation shows that 

Ai = 0.6218... 

By Proposition IV. 61 the periodic points g„ disappear 
as A increases from Ai to A2. The point qi is the last to 
disappear for a value of A close to A2. All points g„ are 
contained in the light-colored trapping region depicted in 
Fig. Ha). 

Proposition V.8. We have W^{qn) fl B ^ for every 
< A < Ai and every n sufficiently large. In particu- 
lar, the homoclinic class of qn is a transitive horseshoe 
provided that n is sufficiently large. 

Proof. Let Qn = (s„,0„) where s„ and 0„ are given by 
([3]). Hence as 77, —> 00 these points approach the hori- 
zontal line 6 = ttX{1 — A)/2. It is also straightforward 



to check that W^{qn) contains the horizontal segment 
joining (0, 6'„) to g„. Let 

00 

a{e) = 1 -^tan(A*6') 

i=0 

be the map whose graph is Soo, the upper bound of B. 
This function is decreasing in A, and by definition of Ai 
we have (T(7rA(l — A)/2) — when A = Ai. Also, we 
must have a(7rA(l — A)/2) > for every < A < Ai. 
Thus, a{dn) > for all large enough n. This proves that 
VF"(g„), which contains the segment joining (0, 9n) to g„, 
intersects the set B C W^{P), bounded from above by 
the graph of cr. Hence this homoclinic class is a horse- 
shoe. To complete the proof, we just need to observe 
that transitivity is a general property of the homoclinic 
classes (see Ref.*? (Ch. IX, Prop. 5.2)). □ 

Proposition V.9. The following statements hold for 
A > Ai.- 

1. W'{P) = BU(j)^\B). 

2. T ^ A/+ \ (P) is a trapping region, i. e. (px (F) C 
int(r). 

3. W{qn) r\B = (D for every n>l. 

Proof. Since the second and third statements immedi- 
ately follow from the first one, we only prove that 

W%P) ^ B[J(j}-^{B) . 

By definition, W'(P) = U„>o'^a"(^)- Thus, 

BU(j)^\B) C W%P). 

To prove the opposite inclusion it is sufficient to prove 
that 

c^f{B)ciBU<l>^\B). 

Suppose that the previous inclusion does not holds, i.e. 
there exists x G (f)^'^{B) such that neither x £ B nor 
X e cf)x\B). Thus (j)x{x) e (f>x\B) \ B. On the other 
hand, it is clear that for every A > Ai we have 

i4>^\B)\B)ncPxiM)^9, 
yielding a contradiction. □ 

Proposition V.IO. There exists 6 > such that the 
periodic points qn are all homoclinically related with px 
for every < X < S. 

Proof. For A > close to 0, the local stable manifold 
^ioc(P^) the graph of the concave monotonic function 
hx connecting the left side s = to the vertex {s,9) = 
(1,0). Using a graph transform argument we can prove 
that the local stable manifold Wf^^{qn) is the graph of a 
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concave monotonic function s — hn{0) with o'n+i(6') < 
hn{0) < cr,i(6'), where 

cr„(0) = 1 - (tan(6i) + tan(Ae') + . . . + tan(A""i6l)) . 

Notice that Mf = { {s,9) : cr„+i(6') < s < cr„(6l) } is the 
region of all points in Mi mapped by /" into the domain 
M2- The graph s = hn{0) also connects the left side 
s = to the vertex {s,d) = (1,0). An easy computation 
shows that 

oo 

1 - ^tan(A'6l„) > 

4 = 1 

and 9n > 7rA/4, the second inequality for A < 1/2. Be- 
cause /i(0,7r/2) = (l,7rA/4) we can deduce from the 
inequality 6n > ttX/A that the local unstable mani- 
fold of Qn is the horizontal segment connecting (O,0„) 
to (1,6'„). Whence = [0' 1] x i^n} intersects 

Wf^^{qm) = gi'aph(ft,m) for any pair of integers n,m > I. 
Since p\ = qo, the proposition is proved. □ 

For future use, we now introduce the new constant 

Ao inf{A > : 3n>l, W''{px)nW%qn) = 0} • 

From Proposition lV.lOl it follows that Xq > 6 > 0, where 
6 is as in Proposition I V. lOl Numerically, we found that 

Ao > 0.6104 . 



VI. BIFURCATION OF THE LIMIT SET 

Recall that fix is the nonwandering set of the map 
In this last section, we formulate a conjecture on 
the decomposition of fix, and discuss the changes in this 
decomposition as the parameter A varies. 



Conjecture VI. 1. For any < A < 1, the nonwander- 
ing set fix 'is a union of three sets: 

nx = PUHxDAx, 

where P is the parabolic attractor introduced in Section 
I VA[ Ax is a hyperbolic transitive attractor, and Hx is 
a horseshoe. Moreover, Hx is either transitive or else 
a (possibly empty) hyperbolic chain of transitive horse- 
shoes. In particular, 

M+ = WiP) U W'{Hx) U WiAx). 

Our next conjecture is justified by the fact that P con- 
sists of periodic points, and the set Ax^JEx is hyperbolic. 

Conjecture VI. 2. The set of periodic points is dense in 
fix. 

The rest of the section is devoted to the justification 
of the previous Conjecture IVI.H and to the analysis of 
the changes in the sets Hx and Ax as A varies. The con- 
clusions based on numerical observations are presented 
as conjectures, whereas the conclusions based on analyt- 
ical arguments are presented as propositions with their 
proofs. We split our discussion into four parts, each cor- 
responding to A taking values inside one of the following 
intervals: (0, Ao), (Ao, Ai), (Ai, A2) and (A2, 1). SeeFig.[71 

A. < A < Ao 

The following conjecture is suggested by numerical 
computations of the invariant manifolds for the points 
qn (see Fig.UJa)). 

Conjecture VI. 3. The manifolds VF"(g„+i) and 
W'(g„) intersect transversally for every < A < Cn+i- 
Moreover, provided that < A < Ai, all qn are mutually 
homoclinically related for sufficiently large n. 
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FIG. 8: Zoom of the phase space with maximal local invariant manifolds of the periodic points, (a) Points g„ 
together with their local stable (dashed green) and unstable (dotted red) manifolds for A = 0.6. The black curves 

represent some iterates of the singular set . 
(b) Points Pn together with their local stable (dashed green) and unstable (dotted red) manifolds for 

A = 0.85. 



In light of Conjecture |VL2l the next conjecture simply 
states that rt\\P is the union of the homoclinic classes 
of the Qn. 

Conjecture VI. 4. Suppose that x <E ft\ \ P is a pe- 
riodic point of Then there are m,n > such that 
W^{x) and W^^Qn) intersect transver sally, and W^{qm) 
and W'^{x) intersect transver sally. 

Proposition VI. 5. // Conjectures I V7. and \ VI^ hold, 
then fix = P U H\ and Hx is a transitive horseshoe for 
< A < Ao. 

Proof. By definition of Aq all q„ are homoclinically re- 
lated for < A < Aq. By Conjectures Ell and |Vl31 
il\ = P U Hx, where Hx denotes the homoclinic class 
of px — (Jo- Then Proposition IV. 81 shows that Hx is a 
transitive horseshoe. □ 



B. Ao < A < Ai 

In this parameter range the set fix \ P splits into 
two or more homoclinic classes dynamically partially or- 
dered. At the bottom of this hierarchy of homoclinic 
classes lies a transitive hyperbolic attractor, and at the 
top a transitive horseshoe whose unstable set intersects 
the basin of attraction of P. We write H ^ H' for 
W'''{H')r\W'{H) 7^0. 

Proposition VI. 6. If Coniectures \VI.0^VI.4\ hold, then 
for Aq < A < Ai, there exists N > I such that fix = 
P U Co U Ci U . . . U Cat and 

1. each C'i is the homoclinic class of some periodic 
point qn. , 

2. Ci n Cj = whenever i ^ j , 



3. Co -< Ci -< . . . ^ Cm ) 

4-. Co is a transitive hyperbolic attractor, 

5. P ^ Cn. 

Proof. Let Co, Ci, . . . , Cat be the homoclinic classes of 
the periodic points By conjecture IVI.4i we have 
flxX P ~ Co U Ci U . . . U Cat. These sets are obviously 
disjoint. Conjecture IVI.3I implies the sets Ci are ordered 
in a finite chain, and we can always display them as in 
item 3. Co is the homoclinic class of the fixed point qo — 
Px, and hence a transitive invariant set. It is attracting 
since it lies at the chain's bottom, and it is hyperbolic 
because of Proposition 1111.31 The set = Ci U . . . U Cat 
is a chain of transitive hyperbolic horseshoes. Finally, 
since Cat is at the chain's top. Proposition IV. 81 implies 
that P -< Cat. □ 

By the definition of Aq, for every Ao < A < Ai there 
is some n > \ such that Ty'(qo) n W''{qn+i) — 0, and, 
in view of Conjecture I VI. 31 this implies there is some 
n >_l such that VF"(g„) n W'{qn+i) = 0. Given n > 1, 
let A„ be the bifurcation point where the homoclinic con- 
nection W'^{qn)C\W''{qn+i) breaks down. The numerical 
value given above for Aq was obtained from the following 
dichotomy: for A < Ao almost every point is attracted to 
P, while for A > Aq there is a non trivial hyperbolic at- 
tractor with an open basin of attraction. We did not try 
to understand these heteroclinic connection breaking bi- 
furcations A„ , but numerical plots indicate that Ao = A„ 
for some rather small n, probably n < 3. 

C. Ai < A < A2 

In this parameter range the periodic points (/„ vanish 
one by one. More precisely, according to Proposition |V|6] 
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there is a decreasing sequence of bifurcation parameters 

Ai < . . . < Cn+l < C„ < . . . < C2 < Ci < A2, 

and qn persists for A < c„, but vanishes for A > c„. 
Hence, unUke the previous interval, only finitely many 
Qn persist for each Ai < A < A2. 

Proposition VI. 7. If Conjectures WI.S^ VI. 4\ hold, then 
for Ai < A < A2, there exists N > such that fix — 
P U Co U Ci U . . . U Ctv and 

1. each Ci is the homoclinic class 0/ some periodic 
point Qn, , 

2. Ci nCj — 9 whenever i ^ j, 

3. Cq ^ C'l -< . . . ^ C'n , 

4- Cq is a transitive hyperbolic attractor, 

5. P T< C, for alli = 0,l,...,N. 

Proof. Keeping the notation of last section, the proof 
here is a simple adaptation of that of Proposition IVL6I 
As before, the hyperbolic attractor Cq is the homoclinic 
class of the fixed point qq = p\. The main differ- 
ence is that for A > Ai, by Proposition IV. 91 we have 
W^{P) = B U 4''^^{B) and there is a trapping region F, 
disjoint from W'^{P), forward invariant under which 
contains all periodic points qn. This proves item 5. □ 

We found numerically that 

• nx\P = Co, for ci < A < A2, 

• flx\P ^ CqUIqi} with qi qo, for C2 < A < ci, 

• nx\P = CoU{(7i}U{(j2} with qi qo and 52 7^ qi, 
for C3 < A < C2. 

D. A2 < A < 1 

By Proposition lV.il in this parameter range the shad- 
owed region A in Fig. [5fb) is a trapping region, i.e. A is 
forward invariant under i/ia- Moreover, all periodic points 
Pn must lie inside A whenever they exist. 

Our numerical analysis suggests the following conjec- 
tures (see Fig. ISKb)). 

Conjecture VI. 8. The periodic points p„ with n < 16 
persist for A2 < A < 1, while those with n > 17 persist 
for A € (0, a„] U [5„, 1), where the sequences bounding the 
gap satisfy a„ \j A2 and fe„ 1. In particular, for any 
given A2 < A < 1, only finitely many points p„ persist. 

Conjecture VI. 9. The periodic points p„ generate two 
homoclinic classes 

1. Co the homoclinic class of the pn with n < 16 or 
A < a„, 

2. Ci the homoclinic class of the Pn with n > 17 and 
A > 6„. 



For \ > 617 (i.e. when Ci becomes non-empty), Cq ^ Ci. 

The next conjecture simply states that fijsj \ P is the 
union of the homoclinic classes of the periodic points pn 
and the fixed point px. 

Conjecture VI. 10. For every A2 < A < 1 and every 
periodic point x G fix \ {P Cl {px}), there exist n,m > 1 
such that W"{x) and W^{pn) intersect transversally, and 
W"{pm) and W^'*(a;) intersect transversally . 

The proof of the following proposition is similar to 
Proposition IVI.6I 

Proposition VI.ll. // Conjecture I VI. 2\ and Conjec- 
tures rVrMVIia hold, then for every A2 < A < 1, we 
have = P U Co U Cl U {px} and 

1. Co is a transitive hyperbolic attractor, 

2. Cl is a transitive horseshoe (possibly empty), 

3. Co n Cl = 0, 

4. Co<Ci< {px}. 
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